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We use a non-detection in v = 1.4 GHz Green Bank Telescope observations of the ultra-faint 
dwarf spheroidal galaxy Segue I, which could be immersed in a non-negligible halo magnetic held 
of the Milky Way, to place bounds on particle dark matter properties. We model the galaxy 
using an Einasto dark matter prohle, and compute the expected synchrotron hux from dark matter 
annihilation as a function of the magnetic held strength B, diffusion coefhcient Do, and particle 
mass for different annihilation channels. The data strongly disfavor annihilations to e^e“ for 
< 50 GeV, but are not sensitive to the bb channel. Adopting a hducial H ~ 2 pG inferred from 
Segue Fs proximity to the Milky Way, our models of annihilation to with = 30 GeV require 
an intermediate value of Dq for consistency with the data. The most compelling limits are obtained 
for WIMP annihilation to p~-. we exclude < 30GeV —>■ p'^p~ at 95% conhdence, unless Dq 
exceeds the Milky Way value or B is signihcantly smaller than we have assumed. 

PACS numbers: 95.35.+d, 98.52.Wz, 98.56.Wm 


I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) are one 
of the leading candidates for the dark matter of the Uni¬ 
verse, and a huge world-wide effort is underway to detect 
them through direct, indirect, and collider experiments. 
WIMPs are natural dark matter candidates because they 
were put forward to solve problems in particle physics 
unrelated to the dark matter puzzle. They also predict 
the correct relic density at the present epoch almost inde¬ 
pendently of the mass. More importantly, the presence of 
weak interactions allows us to make testable predictions. 

Indirect dark matter detection experiments include 
gamma ray observations of the Milky Way center and 
dwarf galaxies mm, cosmic microwave background mea¬ 
surements IM], Cherenkov observations [8] and neutrino 
flux measurements [giiQ]. Interestingly, observations of 
the Milky Way center m by the Fermi gamma ray tele¬ 
scope seem to indicate an excess of gamma rays consis¬ 
tent with dark matter particles of mass = 31 — 40 
GeV annihilating at a rate (fTaw) = (1.4 — 2.0) x 10“^® 
cm^/s, close to the thermal relic value. More recently, 
authors [HUS] found an excess of gamma rays from the 
dwarf galaxy Reticulum II, possibly consistent with dark 
matter annihilation. Models favour very light WIMPs for 
annihilation into and |12j . with most likely 

values of my_ ~ 15 GeV and ^10 GeV for these chan¬ 
nels respectively. 

In this article, we explore dark matter detection 
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through synchrotron radiation from dwarf galaxies. Par¬ 
ticle annihilation results in the release of high energy 
charged particles which emit synchrotron radiation in a 
magnetic field. The detection or absence of this syn¬ 
chrotron radiation can be used to place useful bounds on 
dark matter properties [HI [15] . 

Dwarf spheroidal galaxies are interesting targets for 
synchrotron radiation searches from annihilating dark 
matter because of their proximity and their inferred high 
dark matter content [TBHT^ . In particular, observed line- 
of-sight velocities of individual stars belonging to the 
ultra-faint dwarfs suggest that they have significantly 
higher mass-to-light ratios within their half-radii than 
the classical dwarfs [3D], making them the most likely 
places to observe an electromagnetic signature of anni¬ 
hilating dark matter. Some of the extreme ultra-faint 
dwarfs {L < lO^L©) have a well-measured velocity dis¬ 
persion, and the appearance of being in dynamical equi¬ 
librium [2TH2^ . 

Several radio synchrotron searches for annihilating 
dark matter signatures in dwarf spheroidal galaxies have 
been carried out. In previous work [211 hereafter Pa¬ 
per I], our group obtained deep v = 1.4GHz radio ob¬ 
servations of Draco, Ursa Major II, Goma Berenices, and 
Willman I using the Robert G. Byrd Green Bank Tele¬ 
scope (GBT), and used the resulting maps of Ursa Major 
II and Willman I to constrain models considered by m- 
A more detailed analysis of dark matter annihilation in 
Ursa Major II was described in [211 hereafter Paper II]. 
More recently, authors [2SH2HI used the Australia Tele¬ 
scope Gompact Array (ATCA) to search for radio emis¬ 
sion from the dwarf galaxies visible from the southern 
hemisphere. These observations place tighter constraints 
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FIG. 1: GBT V = 1.4 GHz radio continuum maps in the Segue I region. Plot (a) shows the calibrated GET map, which 
includes background point sources and diffuse sky emission. The linear radiation specific intensity scale ranges from -10 to 
250mJy/beam. Plot (b) shows the mean-subtracted residual map, with point sources removed using the NVSS catalog and 
diffuse sky removed using Planck maps. The linear radiation specific intensity scale ranges from -15 to 30mJy/beam. The 
horizontal line in plot (a) is 1 degree in length, and applies to both panels. The cross in both plots denotes the optical centroid 
of Segue I, and the symbol size resembles the 9 arcmin elliptical half-light diameter of the dwarf. We find no excess radio 
emission at the location of Segue I in the mean-subtracted residual map. 


than in Paper I and Paper II for compact annihilation 
signatures, but are not sensitive to the radio emission on 
scales ^ 15 arcmin produced by many annihilation chan¬ 
nels. 

Synchrotron radiation searches for annihilating dark 
matter hinge on the magnetic field strength B of the 
source, but not much is known about B in dwarf 
spheroidal galaxies [53] ■ The previous work described 
above relied on scaling arguments from detections in 
dwarf irregular galaxies |3()j or optimistic extrapolations 
of the Milky Way’s vertical field m in adopting B ^ 
IfiG. However, our recent search for atomic gas in dwarf 
spheroidal galaxies |32] suggests that they have a more 
tenuous interstellar medium than expected from stellar 
evolution models. This raises the possibility that B is 
much lower than the fiducial value previously invoked, 
producing correspondingly weaker particle dark matter 
constraints. 

Segue I (RA = 10:07:03.2 ± 1.7^ DEC = -bl6:04:25 
± 15”) [3S] is an ultra-faint dwarf spheroidal galaxy at 
a distance ^ 23 kpc from the sun, with a mass-to-light 
ratio M/L ~ SAOO Mq/Lq [5D1 |^, making it one of 
the darkest galaxies known so far. The large mass-to- 
light ratio combined with the close proximity of Segue I 
make it an excellent target for indirect dark matter 
searches. More importantly. Segue I is among the clos¬ 
est dwarf spheroidal galaxies to the sun, which means 
that it could be immersed in the outer halo magnetic 
field of the Milky Way itself [27]. We therefore expect a 
non-negligible ambient magnetic field at the location of 
Segue I, which makes it a unique target for synchrotron 


radiation searches for annihilating dark matter. In this 
article, we make use of new GBT observations of Segue I 
to carry out this search. 


II. GBT OBSERVATIONS AND ANALYSIS OF 
SEGUE I 

Observations sensitive to degree-scale radio emission 
from dwarf spheroidal galaxies require the use of a sin¬ 
gle dish telescope. On the other hand, the low angular 
resolution of single dishes means that we require an in¬ 
terferometer to identify and subtract background point 
sources. The GBT is a fully steerable 100m single dish 
antenna in West Virginia that is well-suited to study 
radio emission in the 300 MHz - 100 GHz frequency 
range. In order to search for extended radio emission 
from Segue I we used the GBT to map a 4° x 4° square 
region in its vicinity at an observation frequency v = 
1.4 GHz, affording the use of publicly available, 45 arc- 
second resolution NRAO VLA Sky Survey (NVSS) data 
[35] at that same frequency to carry out the point source 
subtraction. 

The GBT observations of Segue I were obtained in 
on-the-fly mode during several sessions in late 2013 un¬ 
der program AGBT13B253 using the same setup as de¬ 
scribed in detail in Paper I. The maps are calibrated, 
point-source and diffuse-sky subtracted in a single step by 
building models of the point source and diffuse sky from 
the NVSS catalog and the Planck 30 GHz synchrotron 
model [35]. Both NVSS and Planck templates are con- 
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FIG. 2: Synchrotron radiation from dark matter annihilation. Plots (a) and (b) shows the synchrotron kernel Psynch as a 
function of observation frequency v and magnetic field strength B, respectively. The sharp cut-off at low particle energies 
Pcutoff implies a similarly sharp cutoff at small WIMP masses < Pcutoff. Plot (c) shows the particle spectrum from dark 
matter annihilation for different annihilation channels, and plot (d) shows the particle number density density per unit energy 
after transport at r = 0.1 kpc, for mj,. = 50 GeV, B — 2 fj,G, and Do = 1 kpc^/Gyr for those same annihilation channels. The 
synchrotron radiation at a given observing frequency is sensitive to B, and the annihilation channel. 


volved with the GBT beam and simultaneously fit to the 
GBT map, with a single amplitude for each template 
map. The amplitude of the fit determines the calibration 
of the GBT maps into radiation specific intensity units. 

We use an improved mapping technique to recover 
more extended radio emission (and thus a higher sen¬ 
sitivity to dark matter annihilation signals) in the final 
maps than achieved in Paper I. Single dish data, particu¬ 
larly continuum data, have time variable mean levels due 
to varying ground spillover, changing atmospheric emis¬ 
sion, and low-level electronic gain fluctuations. These 
temporal baselines were removed in Paper I by fitting 
simple linear functions to each scan, which had the side- 
effect of suppressing extended structure in the map. The 
effect is illustrated in fig. 7 of Paper II, where as much 
as 85% of the extended map flux is filtered out. We 
have therefore implemented a de-striping calibration al¬ 
gorithm [37H3S1 that uses only the data differences to 
solve for the baseline terms, resulting in full sensitivity 


to diffuse sky structure. 

Fig.[T] shows the calibrated GBT maps of Segue I in 
(RA, DEG) coordinates in standard radiation specific in¬ 
tensity units of mjy per 9.1 arcmin telescope beam full- 
width at half-maximum diameter (mJy/beam). The hor¬ 
izontal line in plot (a) is 1 degree in length, and applies to 
both panels. The crosses in both panels show the stellar 
centroid location of Segue I, with the symbol size reflect¬ 
ing the ~9 arcmin elliptical half-light diameter of the 
dwarf. Plot (a) shows the calibrated GBT map, which 
includes background point sources and diffuse sky emis¬ 
sion. The latter is most evident in the top-left corner of 
Fig.[TJ where a 1.5-deg “warm spot” clearly underlies the 
background point sources. 

Plot (b) of Fig. shows the mean-subtracted resid¬ 
ual map, now on a much narrower linear intensity scale 
than in plot (a), after background point-source removal 
using the NVSS data and diffuse sky emission removal 
using the Planck map. Some point-like residuals remain 
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FIG. 3: Predicted u — lA GHz radio emission from dark matter annihilation in Segue I for Bin (0° < 6 < 1.12°) for different 
annihilation channels, as a function of the WIMP mass m^. All models adopt a fiducial magnetic field B = 2pG. Plot (a) 
shows models with a diffusion coefficient Do = Ikpc^/Gyr, and plot (b) shows the same models but with Do = lOkpc^/Gyr. 
The dashed horizontal line in both plots shows the RMS fluctuation for Bin ^1 in the mean-subtracted residual map in Fig. 
obtained from simulated maps with the same power spectrum and variance as the data. The bb channel is undetectable in our 
data, whereas our non-detection in Bin ^1 strongly rules out light WIMPs annihilating into e''"e“. 


in the vicinity of the brighter sources; this could indicate 
variability in the ~20-year time span between the NVSS 
and GBT data acquisition, or could arise from small as¬ 
trometric differences between these datasets. These arti¬ 
facts are on a much smaller angular scale than the dark 
matter annihilation signatures we are searching for, and 
are therefore unimportant in our analysis. There are also 
some larger-scale features in the map in Fig. [^, such as 
the cold spot to the West of Segue I and the warm spot 
to the North. They likely represent residual Milky Way 
foreground emission that was not adequately subtracted 
with the scaled Planck map. This is not surprising, since 
small spectral index variations in these features between 
the ;/ = 30 GHz Planck observation frequency and our 
iz = 1.4 GHz observations would lead to improper sub¬ 
traction using our adopted technique. As in Paper I, 
these foregrounds limit the ultimate map sensitivity. 


Fig.^ shows that the optical centroid of Segue I sits 
on the edge of a cold spot in the mean-subtracted residual 
map. In particular, we find no residual emission that 
spatially correlates with Segue I. We therefore proceed to 
compare this non-detection to the synchrotron radiation 
signatures expected for different dark matter annihilation 
channels. 


III. SYNCHROTRON RADIATION FROM 
DARK MATTER PARTICLE ANNIHILATION 


We model the synchrotron radiation from dark mat¬ 
ter particle annihilation in Segue I in a similar manner 
as in Paper H. The transport of charged particles is af¬ 
fected by diffusion and energy losses. Under the sim¬ 
plifying assumptions of spherical symmetry and a con¬ 
stant magnetic field, we may solve for the number of elec¬ 
trons/positrons per unit energy per unit volume dn/dE 
at a distance r from the center of the dwarf [Hill]: 

^ f dn\ d f dn\ „ 


The source term Q is assumed to be entirely due to dark 
matter annihilation: 


Q(r,U) 


{<Ja.v) 2 ^ 

dE 


( 2 ) 


where (aa'c) is the annihilation rate assumed constant at 
the thermal value in all our models |40| : 


(^au) = 2.2 X 10 cm^/s , (3) 


and dN/dE is the number of e“'"e particles produced 
per WIMP: 




( 4 ) 
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TABLE I: Binned data and model radiation specific intensities 


Bin 

^data 

Emodsi, Do = Ikpc’^/Gyr 
rrix = 10 GeV = 30 GeV 

T^T~ P^ p~ r'^T~ p^ P~ 

Dmodei, Do = lOkpc^/Gyr 
= 10 GeV = 30 GeV 

T^r~ p^ p~ T^T~ p^ p~ 

(mjy/arcmin^) 
^2 (mjy/arcmin^) 
#3 (mjy/arcmin^) 
#4 (mjy/arcmin^) 

0.010 ±0.013 

0.018 ±0.017 

-0.0097 ±0.017 

0.0029 ±0.013 

0.026 0.119 0.015 0.037 

0.0025 -0.0089 0.0002 0.00015 

-0.010 -0.049 -0.0061 -0.0152 

-0.013 -0.062 -0.0081 -0.022 

0.0037 0.015 0.0022 0.0046 

0.0037 -0.0017 0.000016 0.00003 

-0.0015 -0.0067 -0.00091 -0.0021 

-0.0024 -0.0064 -0.0015 -0.0023 


Note that we have chosen dN/dE to be the particle 
spectrum per WIMP, and not per annihilation (otherwise 
there would be an extra factor of 2 in the denominator 
of Eq. . The form of dN/dE depends on the annihila¬ 
tion channel. For direct annihilation of WIMPs to 
dN/dE resembles a delta function at the WIMP mass. 
For other channels, dN/dE is a continuous function of 
energy. 

The dark matter density profile of Segue I is as¬ 
sumed to follow the Einasto form: 


Px(^) = P^exp 



(5) 


We choose a = 0.3 to obtain an inner density slope p ^ 
j,-o .3 agreement with The HI Nearby Galaxy Survey 
(THINGS) mass models [1TH15] . We choose Vg = 0.15 
kpc |44l|43 and ps = 6.6 GeV/cm^ in order to reproduce 
the Fermi computation [1] of the emission measure for 
Segue I: 



where the integral over solid angle is from 0 to 0.5 de¬ 
grees, s is the distance measured along the line of sight, 
and the assumed distance to Segue I is L = 23 kpc. 

In Eq. I? is the diffusion parameter which takes the 
form: 


E(E)=Do(^^y , (7) 

where Eq is the diffusion coefficient, and the index 7 = 
0.7 in accordance with the median Milky Way value |3B]. 
b is the energy loss term |14) : 


, ( 8 ) 

where Bq = 0.788[1 -I- 0.102(i3/Ho)^] GeV/Gyr [2] with 
Bq = 1 pG and Eq = 1 GeV. With our assumptions, Eq. 
has a closed form solution [HI [15]. 

Not much is known about the magnetic fields of dwarf 
spheroidal galaxies, although our recent atomic hydrogen 


upper limits [32] suggest that B therein may be negligi¬ 
bly small. However, the magnetic field of the Milky Way 
itself, in which Segue I is bathed by virtue of it proxim¬ 
ity to the sun, is better understood. Authors [47] find 
a magnetic field strength B ^ 6.1 pG at the sun’s loca¬ 
tion consistent with the radio synchrotron background at 
i' = 408 MHz. Authors |3T] used Fermi telescope obser¬ 
vations to detect two giant radio lobes near the Galactic 
bulge, with fields up to 15 pG, extending up to 7 kpc 
from the center. Such a large magnetic field seems in¬ 
consistent with an exponential fall off with distance from 
the center, as commonly assumed [e.g. [48]. Authors [27] 
have suggested a linear scaling. If this is true, then B in 
the outer halo is substantial, and can be non-negligible 
at the location of the dwarf galaxies in the neighborhood 
of the Milky Way. Based on this linear scaling. Segue I 
at L = 23 kpc from the sun is bathed in a magnetic field 
of strength B ^ 2.2 pG. We therefore consider B = 2 pG 
as a fiducial value for Segue I. 

The diffusion coefficient Dq is more difficult to estimate 
for dwarf spheroidal galaxies. Authors iSlEO] adopted 
Dq (X V X I, where v and I are the characteristic veloc¬ 
ity and length scale of the gas of stochastic gas motions. 
From the analysis of peaked iron abundance profiles in 
several galaxy clusters, they found that Dq ~ 10^® cm^/s 
~ 0.3 kpc^/Myr therein [^. This value exceeds the one 
for the Milky Way obtained by analyzing the ratio of 
Boron and Carbon isotopes (T>o,mw = 0.01 kpc^/Myr, 
[46] 1 by over an order of magnitude. Authors m have 
suggested that if this scaling extends to the dwarf galaxy 
regime, then a comparison of the virial velocity disper¬ 
sions of the Milky Way and typical ultra-faint dwarfs 
implies that Dq in the latter systems is an order of mag¬ 
nitude smaller than I?o,mw- We therefore choose our 
fiducial values of Dq to be in the range 1-10 kpc^/Gyr. 

Once Eq. [^has been solved for dn/dE, we may com¬ 
pute the local emissivity: 

r^x dYi 

Jsynch(^5^) — / ^E Tsynch ■ (9) 

J rrie 

The synchrotron kernel Psynch(*^,F) is given by [l4l[52] : 


T^synch — 


2 .70 

GeV 


2.3 


tqcB / d9 sin^ 9 F ( ) 

./n V sm 6» / 

d9siT/^9F (^) , ( 10 ) 

V sin 0 ) 


GHzGyr \fiG J Jq 
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FIG. 4: Exclusion curves as a function of diffusion coefficient Dq and magnetic field strength B for dark matter annihilation to 
l-i~ and T^T~ and WIMP masses rrix = 30GeV in plot (a) and rrix = SOGeV in plot (b). Models with {B, Do) below the 
curves are excluded at 95% confidence. Our models of the channel require only an intermediate value of Do {Do ^ 5.5 

kpc^/Gyr) for our fiducial B ~ 2pG and > SOGeV for consistency with the data. Our most compelling limits come from 
the channel: for = 30 GeV, our non-detection implies that Do ^ 10 kpc^/Gyr - larger than the Milky Way value - 

for B > 1.7 pG at 95% confidence. 


where rg = ^ is the classical electron radius, x is IV. RESULTS 

given by 


X 


0.87 


V 

1.4 GHz 



10 GeV y 


( 11 ) 


and F{t) is computed using [Tlj : 


F{t) = t 

( 12 ) 

The specific intensity X (energy per unit time per unit 
area per unit frequency per unit solid angle) may then 
be calculated by integrating jsynch along the line of sight 
to the dwarf galaxy: 

X(r') = - j ds Jsynch(^; ■ (1^) 

Ji.o.s. 

Fig. shows the form of the synchrotron kernel Psynch 
as a function of particle energy for different observation 
frequencies (plot (a)), and for different values of the mag¬ 
netic field (plot (b)). Psynch das a sharp cutoff at low 
energies which implies a corresponding cutoff at low par¬ 
ticle masses = Pcutoff- Plot (c) shows the energy 
spectrum at the point of annihilation, for different chan¬ 
nels. Direct annihilation to e“'"e“ results in a spectrum 
that is close to a delta function, with most particle at en¬ 
ergies « rriy-. The hadronic channel bb results in a much 
more gradual variation of the spectrum with particle en¬ 
ergy. Plot (d) shows dn/dE after transport, for = 50 
GeV, B = 2 fiG, Dq = Ikpc^/Gyr, at position r = 0.1 
kpc, for the different annihilation channels. 


/ OO 

dzK^fsiz) « 1.25t^/^e-* [648+ 


The mean-subtracted residual radio continuum map 
presented in m shows no excess radio emission at the 
location of Segue I. We now quantify the implications 
of this non-detection for particle dark matter properties 
using the models in §III| To maximize our ability to dis¬ 
tinguish between the dark matter synchrotron emission 
and small-scale noise or artifacts, we divide the map into 
four ~ 1 degree-sized bins centered on Segue I that con¬ 
tain an equal number of pixels: Bin (0° < 9 < 1.12°), 
Bin #2 (1.12° <9 < 1.86°), Bin #3 (1.86° <9 < 2.55°), 
and Bin ^4 (2.55° < 9 < 4.0°), where 6* = 0° corre¬ 
sponds to the location of Segue I. In order to compute 
the root mean square (RMS) fluctuation in intensity per 
bin, we use the LensTools package [S3] to simulate 1000 
maps with the same power spectrum and variance as our 
observed map, and compute the mean intensity and RMS 
in each bin. 

Fig. i shows the radiation intensity in Bin ^1 (closest 
to the Segue I stellar centroid) as a function of the WIMP 
mass with B = 2 /iG, for Dq = 1 kpc^/Gyr in plot (a) and 
Dq = lOkpc^/Gyr in plot (b). Also shown with dashed 
lines is the RMS fluctuation in the data Bin that 
we obtained from the simulations. We note that the bb 
channel is undetectable, since it is well below the noise 
level even for the optimistic choice Dq = Ikpc^/Gyr. 
The leptonic channels are generally easier to constrain, 
particularly for small Dq and low m^. The e“'"e“ channel 
predicts the largest Bin ^1 intensity for small WIMP 
masses, but falls off very quickly with rising mass. This 
is because the energy spectrum of e+e” due to WIMP 
annihilation is close to a delta function at the WIMP 
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mass (see Fig. [!)• Fig. I shows that we rule out light 
WIMPs that annihilate to e+e” for Dq < lOkpc^/Gyr, 
a model that is also strongly disfavored by the polarized 
cosmic microwave background [7] . We therefore focus on 
the constraints obtained for the and fJ.~ channels 
below. 

The radiation intensity values in the four bins due to 
dark matter annihilation are computed for the and 

annihilation channels, and for different choices of 
B, Dq, and m^. As shown by Fig. the synchrotron 
emission from dark matter annihilation in Bin will 
be large and positive. Bin ^4 is expected to be negative 
since the residual map in Fig. is mean-subtracted. 
Bins #2 and #3 may be either positive or negative de¬ 
pending on the annihilation channel, values of Dq, etc. 
On the other hand, in the absence of a detection we do 
not expect the data to show any definite trend in inten¬ 
sity values as we move from Bin to Bin ^4. 

Table 1 gives the observed mean-subtracted residual 
map intensity values and error bars in Bins #1, #2, ^3 
and ^4, along with theoretical predictions for WIMP 
mass = 10 and 30GeV and for Dq = 1 and 10 
kpc^/Gyr for the and channels. The mean 

intensity of the data in all four bins is consistent with 
zero within uncertainty. As expected from Fig.[^ Bin ^1 
values for the models are positive and larger for smaller 
WIMP masses and smaller values of Dq. Model Bin 
values are all negative. By and large, the model values 
in Table 1 exceed those for the data in (absolute) am¬ 
plitude in all bins, particularly for = lOGeV; we 
can therefore place competitive bounds on the r+r“ and 
annihilation channel WIMP properties for Dq < 10 
kpcVGyr and S < 2/rG. 

We now determine the 95% exclusion curves in the 
Dq — B plane for WIMPs of different masses annihilating 
into r+T” and by computing the quantity for 

each model, summed over the four bins: 


X 


2 



(14) 


where ai is the RMS fluctuation in each bin obtained 
from the data simulations. The likelihood function is 
then given by: 

/: = exp|-i(x2-X%„)| . (15) 

Each exclusion curve is obtained by computing the cen¬ 
tral 95% area under the likelihood curve, assuming uni¬ 
form priors. 

Fig. 1^ shows 95% exclusion curves for WIMP annihila¬ 
tions into r+T“ and for two representative WIMP 

masses, = 30GeV in plot (a) and = 50GeV 
in plot (b). The region below each curve is excluded 
by the data. Recall that the fiducial values for the 
diffusion coefficient and the magnetic field strength are 
Ho = 1 — 10 kpc^/Gyr (the lower bound from scaling ar¬ 
guments, and the upper bound equal to the Milky Way 


value) and B = 2 fiG (given Segue Ps proximity to the 
Milky Way), respectively (®- The exclusion limits are 
therefore quite strong for light WIMPs annihilating to 
for = 30 GeV, our non-detection implies that 
Dq > 10 kpc^/Gyr - larger than the Milky Way value - 
for B > 1.7 fiG, while we exclude = 30 GeV —?> 
for i? > 1 ^G and Dq ^ 2 kpc^/Gyr. In combination 
with the values of Table 1, Fig. [^illustrates that WIMPs 
with ^ 30 GeV annihilating to - such as those 

favored for that channel in the dark matter interpreta¬ 
tion of the Reticulum II gamma-ray excess m - are ruled 
out by our data unless B is significantly smaller (or Dq 
is significantly larger) than argued here. 

Fig. i shows that our limits on the channel 

are weaker, requiring only that Dq > 5.5 kpc^/Gyr for 
H ~ 2 fiG and > 30 GeV for consistency with the 
data. It is therefore unlikely that we would detect the 
= 30 GeV —>■ r+r” signature if Dq ~ 10 kpc^/Gyr, 
as might be expected if the synchrotron radiation in 
Segue I is generated by the magnetic field of the Milky 
Way in which it is bathed. 


V. CONCLUSIONS 

In this article, we placed competitive bounds on dark 
matter properties using radio observations of the ultra¬ 
faint dwarf spheroidal galaxy Segue I. This target is par¬ 
ticularly compelling for radio synchrotron dark matter 
searches, because its close proximity (L ^ 23 kpc) implies 
that it could be bathed in a non-negligible halo magnetic 
field of the Milky Way. (jn] presented a ly = 1.4 GHz map 
of a 4° X 4° region around Segue I obtained with the 
GBT. Point sources were subtracted from the map using 
the NVSS catalog, and the diffuse sky was subtracted 
using a Planck map. Our non-detection near the stel¬ 
lar centroid of Segue I in the resulting mean-subtracted 
residual map was used to place bounds on the WIMP 
mass for different values of magnetic field B and the dif¬ 
fusion coefficient Dq. 

m presented our models for the radio synchrotron 
emission expected from Segue I for different dark mat¬ 
ter annihilation channels, adopting an Einasto density 
profile for the dwarf and matching its Fermi emission 
measure. We studied the diffusion of charged particles 
in a magnetic field, computing the annihilation rate after 
transport determined by Dq. The synchrotron radiation 
at a given observing frequency is sensitive to H, and 
the annihilation channel. 

§IV|used these models to place bounds on particle dark 
matter properties from our observations. Since the dark 
matter synchrotron intensity has very little power on 
small scales, we binned our data into four ^ 1 degree¬ 
sized bins in order to maximize the dark matter signal 
while minimizing contributions from thermal noise and 
small-scale artifacts. We computed the likelihood func¬ 
tion for different WIMP masses, and placed bounds on 
dark matter properties. Our data strongly rule out light 



WIMPs annihilating to e+e”, but are not sensitive to 
the bb channel. Our models of the t'^t~ channel require 
only an intermediate value of Dq {Dq > 5.5 kpc^/Gyr) for 
our fiducial B ~ 2 fj,G and > 30 GeV for consistency 
with the data. Our most compelling limits come from 
the channel: for = 30 GeV, our non-detection 

implies that Dq > 10 kpc^/Gyr - larger than the Milky 
Way value - for i? > 1.7 fiG at 95% confidence. This 
limit is inconsistent with the ^ 10 GeV —> in¬ 

terpretation of the Reticulum II gamma-ray excess unless 
B is significantly smaller than the fiducial value adopted 
here. 

Segue I is a unique target for radio synchrotron dark 
matter searches because of its close proximity to the 
Milky Way, but a joint radio analysis of several promis¬ 
ing dwarf spheroidal galaxies could also place powerful 
constraints on particle dark matter. We have obtained 
~300 hours of wide-held GET mapping observations in 


the vicinity of 6 such systems in order to do just that. A 
detailed analysis of this data is left to future work. 


Acknowledgments 

We thank Beth Willman, Judith Irwin, and Lawrence 
Widrow for helpful suggestions, and acknowledge assis¬ 
tance in the data reduction from Bryan Fichera and 
Melissa Diamond. AN is grateful for hnancial support 
from Queen’s University, and acknowledges hospitality 
of the University of Pennsylvania. KS acknowledges sup¬ 
port from the Natural Sciences and Engineering Research 
Council of Canada. The National Radio Astronomy Ob¬ 
servatory is a facility of the National Science Founda¬ 
tion operated under cooperative agreement by Associated 
Universities, Inc. 


[1] M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, 
L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bech- 
tol, R. Bellazzini, et al.. Physical Review Letters 107, 
241302 (2011), 1108.3546. 

[2] Fermi-LAT Collaboration, ArXiv e-prints (2015), 
1503.02641. 

[3] X. Chen and M. Kamionkowski, Phys. Rev. D 70, 043502 
(2004), astro-ph/0310473. 

[4] N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. D 
72, 023508 (2005), astro-ph/0503486. 

[5] S. Galli, F. locco, G. Bertone, and A. Melchiorri, Phys. 
Rev. D 80, 023505 (2009), 0905.0003. 

[6] S. Galli, F. locco, G. Bertone, and A. Melchiorri, Phys. 
Rev. D 84, 027302 (2011), 1106.1528. 

[7] A. Natarajan, Phys. Rev. D 85, 083517 (2012), 
1201.3939. 

[8] M. A. Sanchez-Conde, M. Cannoni, F. Zandanel, M. E. 
Gomez, and F. Prada, Journal of Cosmology and As- 
troparticle Physics 12, 011 (2011), 1104.3530. 

[9] IceCube collaboration, M. G. Aartsen, R. Abbasi, Y. Ab- 
dou, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, 
D. Altmann, K. Andeen, et al., ArXiv e-prints (2012), 
1212.4097. 

[10] The Snper-Kamiokande Collaboration, :, K. Choi, 
K. Abe, Y. Haga, Y. Hayato, K. lyogi, J. Kameda, 
Y. Kishimoto, M. Miura, et al., ArXiv e-prints (2015), 
1503.04858. 

[11] T. Daylan, D. P. Finkbeiner, D. Hooper, T. Linden, 
S. K. N. Portillo, N. L. Rodd, and T. R. Slatyer, ArXiv 
e-prints (2014), 1402.6703. 

[12] A. Geringer-Sameth, M. G. Walker, S. M. Koushiappas, 
S. E. Koposov, V. Belokurov, G. Torrealba, and N. Wyn 
Evans, ArXiv e-prints (2015), 1503.02320. 

[13] D. Hooper and T. Linden, ArXiv e-prints (2015), 
1503.06209. 

[14] S. Colafrancesco, S. Profnmo, and P. Ullio, Astron. As- 
trophys. 455, 21 (2006), arXiv:astro-ph/0507575. 

[15] S. Colafrancesco, S. Profumo, and P. Ullio, Phys. Rev. D 
75, 023513 (2007), arXiv:astro-ph/0607073. 

[16] M. L. Mateo, Annnal Review of Astronomy and Astro¬ 


physics 36, 435 (1998), arXiv:astro-ph/9810070. 

[17] L. E. Strigari, S. M. Koushiappas, J. S. Bullock, and 
M. Kaplinghat, Phys. Rev. D 75, 083526 (2007), 
arXiv:astro-ph/0611925. 

[18] L. E. Strigari, S. M. Koushiappas, J. S. Bullock, 
M. Kaplinghat, J. D. Simon, M. Geha, and B. Willman, 
Astrophys. J. 678, 614 (2008), 0709.1510. 

[19] L. E. Strigari, J. S. Bullock, M. Kaplinghat, J. D. Si¬ 
mon, M. Geha, B. Willman, and M. G. Walker, Natnre 
(London) 454, 1096 (2008), 0808.3772. 

[20] M. Geha, B. Willman, J. D. Simon, L. E. Strigari, E. N. 
Kirby, D. R. Law, and J. Strader, Astrophys. J. 692, 
1464 (2009), 0809.2781. 

[21] J. D. Simon and M. Geha, Astrophys. J. 670, 313 (2007), 
0706.0516. 

[22] E. N. Kirby, M. Boylan-Kolchin, J. G. Cohen, M. Geha, 
J. S. Bullock, and M. Kaplinghat, Astrophys. J. 770, 16 
(2013), 1304.6080. 

[23] M. Walker, Dark Matter in the Galactic Dwarf Spheroidal 
Satellites (2013), p. 1039. 

[24] K. Spekkens, B. S. Mason, J. E. Aguirre, and B. Nhan, 
ArXiv e-prints (2013), 1301.5306. 

[25] A. Natarajan, J. B. Peterson, T. C. Voytek, K. Spekkens, 
B. Mason, J. Agnirre, and B. Willman, Phys. Rev. D 88, 
083535 (2013), 1308.4979. 

[26] M. Regis, L. Richter, S. Colafrancesco, M. Massardi, 
W. J. G. de Blok, S. Profumo, and N. Orford, Mon. Not. 
R. Astron. Soc. 448, 3731 (2015), 1407.5479. 

[27] M. Regis, L. Richter, S. Colafrancesco, S. Profumo, 
W. J. G. de Blok, and M. Massardi, Mon. Not. R. Astron. 
Soc. 448, 3747 (2015), 1407.5482. 

[28] M. Regis, S. Colafrancesco, S. Profumo, W. J. G. de Blok, 
M. Massardi, and L. Richter, Journal of Cosmology and 
Astroparticle Physics 10, 016 (2014), 1407.4948. 

[29] R. Beck and R. Wielebinski, Magnetic Fields in Galaxies 
(2013), p. 641. 

[30] K. T. Chyzy, M. Wezgowiec, R. Beck, and D. J. Bomans, 
Astron. Astrophys. 529, A94 (2011), 1101.4647. 

[31] E. Carretti, R. M. Crocker, L. Staveley-Smith, 
M. Haverkorn, C. Purcell, B. M. Gaensler, G. Bernard!, 



9 


M. J. Kesteven, and S. Poppi, Nature (London) 493, 66 
(2013), 1301.0512. 

[32] K. Spekkens, N. Urbancic, B. S. Mason, B. Willman, 
and J. E. Aguirre, Astroph. J. Letters 795, L5 (2014), 
1410.0028. 

[33] M. Geha, B. Willman, J. D. Simon, L. E. Strigari, E. N. 
Kirby, D. R. Law, and J. Strader, Astrophys. J. 692, 
1464 (2009), 0809.2781. 

[34] G. D. Martinez, Q. E. Minor, J. Bullock, M. Kaplinghat, 
J. D. Simon, and M. Geha, Astrophys. J. 738, 55 (2011), 
1008.4585. 

[35] J. J. Gondon, W. D. Gotton, E. W. Greisen, Q. F. Yin, 
R. A. Perley, G. B. Taylor, and J. J. Broderick, Astron. 
,J. 115, 1693 (1998). 

[36] Planck Collaboration, R. Adam, P. A. R. Ade, 

N. Aghanim, M. 1. R. Alves, M. Arnaud, M. Ashdown, 
J. Aumont, C. Baccigalupi, A. J. Banday, et ah, ArXiv 
e-prints (2015), 1502.01588. 

[37] C. G. T. Haslam, M. J. S. Quigley, and C. J. Salter, Mon. 
Not. R. Astron. Soc. 147, 405 (1970). 

[38] C. G. T. Haslam, U. Klein, C. J. Salter, H. Stoffel, W. E. 
Wilson, M. N. Cleary, D. J. Cooke, and P. Thomasson, 
Astron. Astrophys. 100 , 209 (1981). 

[39] B. Winkel, L. Elder, and A. Kraus, Astron. Astrophys. 
547, A119 (2012), 1210.4005. 

[40] T. Han, Z. Liu, and A. Natarajan, Journal of High En¬ 
ergy Physics 11, 8 (2013), 1303.3040. 

[41] F. Walter, E. Brinks, W. J. G. de Blok, F. Bigiel, R. C. 
Kennicutt, Jr., M. D. Thornley, and A. Leroy, Astron. 
J. 136, 2563 (2008), 0810.2125. 

[42] W. J. G. de Blok, F. Walter, E. Brinks, C. Trachternach, 


S.-H. Oh, and R. C. Kennicutt, Jr., Astron. J. 136, 2648 
(2008), 0810.2100. 

[43] S.-H. Oh, W. J. G. de Blok, E. Brinks, F. Walter, 
and R. C. Kennicutt, Jr., Astron. J. 141, 193 (2011), 
1011.0899. 

[44] J. Aleksic, S. Ansoldi, L. A. Antonelli, P. Antoranz, 

A. Babic, P. Bangale, U. Barres de Almeida, J. A. Bar¬ 
rio, J. Becerra Gonzalez, W. Bednarek, et al., Journal 
of Cosmology and Astroparticle Physics 2, 008 (2014), 
1312.1535. 

[45] R. Essig, N. Sehgal, and L. E. Strigari, Phys. Rev. D 80, 
023506 (2009), 0902.4750. 

[46] F. Donato, N. Fornengo, D. Maurin, P. Salati, and 
R. Taillet, Phys. Rev. D 69, 063501 (2004), arXiv:astro- 
ph/0306207. 

[47] A. W. Strong, 1. V. Moskalenko, and O. Reimer, Astro¬ 
phys. J. 537, 763 (2000), astro-ph/9811296. 

[48] R. Beck, Astrophysics and Space Sciences Transactions 
5, 43 (2009). 

[49] P. Rebusco, E. Churazov, H. Bdhringer, and W. Forman, 

Mon. Not. R. Astron. Soc. 359, 1041 (2005), arXiv:astro- 
ph/0501141. 

[50] P. Rebusco, E. Churazov, H. Bdhringer, and W. Forman, 

Mon. Not. R. Astron. Soc. 372, 1840 (2006), arXiv:astro- 
ph/0608491. 

[51] T. E. Jeltema and S. Profumo, Astrophys. J. 686, 1045 
(2008), 0805.1054. 

[52] M. S. Longair, High Energy Astrophysics (2011). 

[53] Petri, A., http://www.columbia.edu/~ap3020/LensTools/html 
(2014). 



